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ABSTRACT 

We have obtained time series observations of the Orion Nebula Cluster at 70 /xm 
and 160 ^m from the Herschel /PACS Photometer. This represents the first wide- 
field far-infrared photometric monitoring of a young star forming region. The acquired 
35'x35' maps show complex extended structures, with unprecedented details, that trace 
the interaction between the molecular gas and the young hot stars. We detect 43 proto- 
stars, most of which are situated along the integral-shaped filament extending from the 
Orion nebula, through OMC 2 and to OMC 3. We present high-reliability light curves 
for some of these objects using the first six epochs of our observing program spread 
over 6 weeks. We find amplitude variations in excess of 20% for a fraction of the de- 
tected protostars over periods as short as a few weeks. This is inconsistent with the 
dynamical time-scales of cool far-IR emitting material that orbits at hundreds of AU 
from the protostar, and it suggests that the mechanism(s) responsible for the observed 
variability originates from the inner region of the protostars, likely driven by variable 
mass accretion. 



Subject headings: Stars: formation — Stars: protostars — Stars: Variables: general - 
Infrared: stars 
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Introduction 



Photometric variability was one of the original, defining characteristics of the class of ob- 
jects tha t were lat er determined to be stars in the process of formation, or Young Stellar Objects 
(YSOs) (|Jovlll945l ). While optical and near-IR variability studies have demonst rated the relation 
between accretion sh ocks and hotspots on the surface of rotating YSOs (e.g., IVrba et al.l 119861 ; 



Carpenter et al.ll200ll ). observations at longer wavelengths offer another perspective to study cooler 



circumstellar disks and envelopes surrounding nascent stars. In particular, re cent mid-IR photo 



metric and spectroscopic time series of YSOs with the Spitzer Space Telescope (jWerner et al 



2004) 



have begun to s hed light on inner disk structures, identifying both inner disk warps and 'clouds' 
in the disk (e.g., iMorales-Calderon et aDbood ; iMuzerolle et aDbood ; lEspaillat et all boill). Elab- 
orate models offer satisfying explanations to the observ ed time series ([Flaherty &: Muzerolldl2010l ; 
Turner et allboid : iFlahertv et aLlboill : ke et alJboia l. 



Similar photometric variability studies are a lot sparser in the literature at longer wave l engths 
(A > 50 /im), mostly due to the difficult access to this wavelength domain. lHarvey et al.1 (119981 ) 
report on a 200% flux variation of SSV 13 in NGC 1333 oyer a t wo-year period using the Kuiper 
Airborne Observatory. lJuhasz et al.l (|2007l ) and lSitko et al.l (|2008l ) analyze the 1-300 /im variability 
of the star SV Cep over two years, and they invoke a growing warp at the inner edge of the disk, 
which leads to the shadowing of the outer disk, to reproduce t he observations ac r oss the entire 



wavele ngth range. Based on reprocessed ISOPHOT observations (jLemke et al,lll996l ). iKospal et al 



(|2007l ) present the far-IR fading of the pre-main sequence star OO Serpentis after it went into 
outburst, and they show that the 60 //m flux changed by a factor of 8 over a period of 11 years. 
Most sources targeted by these monitoring programs are disk-bearing YSOs. 

In this article, we present the first sensitive wide-field far-IR photometric monitoring of pro- 
tostars which significantly extends the statistics on YSO variability monitoring in this wavelength 
regime and evolutionary stage. We recently completed the firs t set of time series o bservations of 
the Orion Nebula Cluster (O NC) with the PACS instrument (jPoglitsch et al.ll2010l ) onboard the 
Herschel Space Observatory (jPilbratt et al.ll2010l ). During the first two- month visibility window, 
we obtained six observations of a single field centered on the integral-shaped filament located north 
of the ONC , which contains h undreds of YSOs. In addition, this field has been monitored in 



the optical (IHerbst et al 



(Morales-Calderon et al 



2011 



2000), and in the mid-IR with Spitzer as part of the YSOVAR program 



thus making it the largest YSO variability database to date. 



We first describe the observations and data reduction in section [2j then we give a detailed 
account of the photometric measurements and present the light-curves extracted from the time 
series observations. Finally we discuss the possible origin of the observed variability in section [3l 
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Method 



2.1. Observations 



We use the Herschel/PACS Photometer simultaneously in the blue and red channels, oper- 
ating at 70 and 160 /xm respectively, to measure photometric time series of YSOs in the ONC. 
This Herschel OT1 observing program is designed to make the best use of the observatory while 
accommodating the constraints dictated by the need for precise photometric monitoring. 

The building block of our program is a 30-minute single-direction scan map observation cover- 
ing a 35'x35' field. This scheme of scanning a relatively large area that contains several YSOs was 
preferred over repeated observations of individual sources because of the prohibitive telescope over- 
heads that this would entail. The targeted field was ch osen to contain as man y YSOVAR- monitored 
protostars as possible, with over 100 known in total (jKrvukova et al.ll2012l ). This building block 
is then repeatedly scheduled with absolute timing constraints to cover several 2-months visibility 
windows spread over a period of nearly 2 years. Moreover, within each PACS-Photometer observing 
campaign, which lasts 2.5-days, our observations are scheduled at least 12 hours after the end of 
the cryo-cooler recycling and switch-on of the detectors. This is to ensure the highest thermal and 
electrical stability possible required for photometric variability monitoring. 

In addition, considering that map-making algorithms rely on scan-directio n redundancy to 



filter out instrumental dr ifts from the extended emission present in the field (e.g., ICantalupo et al 



2010 



Roussel et al.ll2012l ). we require that the scan direction of consecutive observations are rotated 



by 90° with respect to the previous one. The scan-direction redundancy is further increased by the 
secular rotation of the spacecraft roll angle along its orbit (~10° between consecutive observations). 
This observing strategy thus allows a very good rejection of instrumental drifts when combining 
single-epoch observations, and consequently the reconstruction of artifact-free extended emission 
structures^ (see Figure [1]). 



2.2. Data Processing 

We retrieve the raw data from the Herschel Science Archive and process individual observations 
to obtain calibrated data cubes (Level 1 frames) using the PA CS standar d pipeline with the Herschel 
Interactive Processing Environment (HIPE version 7.0 1931. lOtt Il2010l ). We correct the bolo meter 
signals for non-linearities based on pre-launch calibration measurements (jBillot et al.ll2006l ). We 
also modify the frame's astrometry to correct for telescope pointing errors (up to 3" in the sixth 
visit). 



1 A Spitzer - Herschel composite image based on this dataset has also been released on the NASA web site at 
http://www.nasa.gov/mission_pages/herschel/multimedia/ pial3959.html 
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We create maps from individual single-direction scan observations using the standard approach 
of highpass filtering the bolometer timelines and projecting the data cube on regular grids with 
pixel sizes of 1" and 2" in the blue and red PACS channels, respectively. We use a filter width twice 
as long as a scan leg (2000 readouts), and we mask bright regions prior to applying the filter to 
minimize shadowing artifacts along the scan direction. This filtering process alters to some extent 
the largest spatial scales in the map but preserves the signal at the scale of point sources. The 
source photometry is measured from these single-epoch maps (see section [2~3l) . 

In parallel, we export individual L evel 1 frames out of HIPE to be combined by the IDL-based 



map-making algorithm Scanamorphos (jRoussel et al.ll2012j ). Figure [T] presents the maps resulting 



from the combination of the six epochs observed during the first visibility window between February 
26, 2011 and April 10, 2011. Most of the YSOs in this field are situated along the Orion A molecular 
clouds OMC-1, -2 and -3, north of the ONC. A small area of 10" radius is saturated in both 
maps at the location of the Trapezium cluster where the signal reaches ~ 4 x 10 6 MJy-sr^ 1 and 
~2x 10 5 MJy.sr -1 in the blue and red channels, respectively. The Orion bar is visible to the south 
east of the Trapezium cluster. A detailed discussion of the extended emission structure in this field 
is beyond the scope of this letter and will be treated in a subsequent article. 



2.3. Source photometry 

The extended emission in the ONC is highly structured (cf Figure [I]) , and it exhibits a great 
dynamic range, in excess of 1:10000. This makes the fine tuning of source-finding algorithms 
very difficult to converge to a satisfying solution. We therefore opt for a more pragmatic visual 
inspection of the combined 70 /xm map to identify point-like sources. We find 43 objects, with fluxes 
ranging from 0.4 to 450 Jy. About 80% of these sources have 160 /xm counterparts. One third of 
those appear to be point-like objects at 160 /xm, with fluxes ranging from 13 to 280 Jy, while the 
other two thirds are spatially resolved and often indistinguishable from the underlying extended 
emission structure. Most of the sources d e tected in the PACS maps have been previously identified 



as Class 1/0 protostars bv lMegeath et al.l (|2012j ). as presented in Table [TJ 



For each single-epoch map, we measure the photometry of point-like sources within a small 
aperture centered on the PSF core (4" and 8" radius in the blue and red bands, respectively). We 
measure background levels, and associated background standard deviations, within an annulus of 
inner/outer radii of 6"/8" and 10"/16" in the blue and red bands, respectively. Although the signal 
in the PSF wings affects to some extent the determination of the background level and standard de- 
viation, it is best to work with small aperture/annulus radii to minimize the effect of the structured 
background emission. We correct the measured fluxes with aperture correction factors derived from 
the PACS PSF@. The PACS Photometer bolometer arrays are extremely stable detectors with gain 



2 The PSFs are measured on the Asteroid Vesta, and are publicly available through the Herschel Science Center 
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variations around 1% over a year timescale ([Billot et al.ll2010l ). Flux uncertainties are thus mainly 



dominated by errors in background level estimates, which are likely underestimated by the standard 
deviation measured in the thin annulus. 

We estimate empirical flux uncertainties as the typical flux variations measured on a given 
source in various combinations of aperture/annulus radii. At 70 /im, we find that flux uncertainties 
are generally below 5%, but they can reach 20% for sources located in regions of highly structured 
backgrounds. At 160 /im, the absolute flux uncertainties are generally higher (~20-50%), mainly 
due to the structured background and to the fact that some sources are extended at this wavelength 
in which case aperture corrections are no longer appropriate and may introduce systematic errors. 
Table [1] gives the source fluxes averaged over the visibility window period. The errors quoted in the 
table reflect the standard deviation measured in the background annulus and likely underestimate 
the flux uncertainties as described above. 

Similarly, we estimate the relative flux uncertainties, i.e. the flux variations measured on 
a given source at different epochs, by exploring various combinations of aperture/annulus radii. 
Presently, if any point of a light curve shows a flux variation of more than 5% in the various 
combinations of aperture/annulus radii explored, the light curve is deemed inaccurate. This rather 
stringent criterion ensures the highest reliability of the light curves. In the sample of 43 sources, we 
find that 17 have a reliable light curve in the blue band, and only 6 in the red band. Table [2] gives 
the time series fluxes of sources that pass this selection criterion, and figure [2] presents a sample of 
these high-reliability light curves. 



3. Results and Discussion 

It is a difficult task to disentangle intrinsic source variability and photometric errors based 
solely on 6 data points, especially for faint sources. Nevertheless, considering the criterions we used 
in selecting reliable light curves (cf section [273]) . we expect that flux variations greater than ~10% 
can be confidently attributed to the intrinsic variability of the observed sources. 

We find that 8 sources out of 17 show peak-to-peak flux variations higher than ~10%. The left 
column of figure [2] presents a sample of such light curves. The case of HOY J053522. 27-050116. 8 
(top left plot) is of particular interest with its sine-like light curve that shows an amplitude of nearly 
20%. The smooth shape of the light curve suggests that the photometric uncertainties are very 
small for this source. In addition both 70 fj,m and 160 /mi light curves show the same trend, which 
strengthens further our confidence that this source is indeed variable in the Far-IR. The object 
HOY J053524. 23-050831. 9 also shows strong variations in excess of 20% in just over 2 weeks. For 
compariso n, we present in figure El Spitzer light curves obtained in 2009 as part of the YSOVAR 



program (jMorales-Calderon et al 



201 ll ). and it appears that objects variable in the far-IR tend 



web pages at http://herschel.esac.esa.int/. 



Table 1. Source flux averaged over the visibility window period. 



Source ID R.A. Dec. Flux (Jy) YSO 

(deg.) (deg.) 160 /im a Class b 



HOY 


J053529 


41 


045851 


6 


83 


872 


-4 


981 


6.33±0.07 


27.4±0.69 


1/0 


HOY 


J053533 


95 


045955 


1 


83 


891 


-4 


998 


0.58±0.01 




1/0 


HOY 


J053534 


22 


050054 


2 


83 


892 


-5 


015 


2.13±0.02 




1/0 


HOY 


J053530 


33 


045938 


6 


83 


876 


-1 


994 


0.93±0.01 




11 


HOY 


J053531 


36 


050015 


9 


83 


880 


-5 


004 


0.48±0.01 




11 


HOY 


J053529 


50 


045952 


7 


83 


872 


-4 


997 


12.9±0.10 


21.6±0.69 


1/0 


HOY 


J053515 


98 


050004 


7 


83 


816 


-5 


001 


7.29±0.08 




1/0 


HOY 


J053518 


12 


050035 


6 


83 


825 


-5 


009 


27.4±0.22 


38.7±1.60 


1/0 


HOY 


J053518 


73 


050053 


1 


83 


828 


-5 


014 


2.82±0.03 




1/0 


HOY 


J053519 


74 


050104 


7 


83 


832 


-5 


017 


2.28±0.05 




1/0 


HOY 


J053522 


27 


050116 


8 


83 


842 


-5 


021 


33.5±0.28 


70.4±1.63 


1/0 


HOY 


J053523 


26 


050132 


6 


83 


846 


-5 


025 


52.5±0.38 


158±2.58 


1/0 


HOY 


J053528 


00 


050343 


1 


83 


866 


-5 


061 


18.0±0.21 


27.2±1.39 


1/0 


HOY 


J053526 


13 


050356 


5 


83 


860 


-5 


065 


72.1±0.75 


73.1±1.61 


1/0 


HOY 


J053519 


61 


050453 


7 


83 


831 


-5 


081 


3.60±0.06 




1/0 


HOY 


J053531 


17 


050549 


1 


83 


881 


-5 


096 


5.63±0.06 




1/0 


HOY 


J053527 


93 


050459 





83 


866 


-5 


083 


2.40±0.02 




1/0 


HOY 


J053525 


89 


050544 


2 


83 


857 


-5 


095 


12.3±0.11 


41.6±1.02 


1/0 


HOY 


J053526 


58 


050610 


7 


83 


860 


-5 


102 


6.80±0.05 




1/0 


HOY 


J053525 


51 


050759 


1 


83 


856 


-5 


133 


10.5±0.09 




1/0 


HOY 


J053527 


56 


050704 


6 


83 


864 


-5 


117 


1.39±0.11 




1/0 


HOY 


J053523 


95 


050753 


7 


83 


849 


-5 


131 


2.88±0.07 




1/0 


HOY 


J053524 


23 


050831 


9 


83 


850 


-5 


142 


7.23±0.10 


18.2±1.33 


1/0 


HOY 


J053524 


71 


051030 


6 


83 


852 


-5 


175 


95. Oil. 21 


36.6±4.19 


1/0 


HOY 


J053526 


91 


051017 


5 


83 


862 


-5 


171 


17.5±0.79 




1/0 


HOY 


J053527 


56 


050934 


5 


83 


864 


_5 


159 


450±6.51 


279±9.42 


1/0 


HOY 


J053523 


29 


051203 





83 


847 


-5 


200 


48.4±0.57 


54.7±2.16 


1/0 


HOY 


J053520 


15 


051317 


4 


83 


833 


-5 


221 


46.8±0.38 




1/0 


HOY 


J053521 


34 


051319 


5 


83 


838 


-5 


222 


6.26±0.09 






HOY 


J053518 


50 


051340 


2 


83 


827 


-5 


227 


3.68±0.12 




11 


HOY 


J053519 


17 


051534 


7 


83 


831 


-5 


259 


27.1±0.59 




1/0 


HOY 


J053517 


93 


051614 


8 


83 


824 


-5 


270 


3.51±0.13 




11 


HOY 


J053531 


42 


052515 


7 


83 


880 


-5 


421 


31.9±0.98 




11 


HOY 


J053515 


69 


052040 


6 


83 


815 


-5 


344 


189±6.60 




11 


HOY 


J053514 


71 


052135 


9 


83 


811 


-5 


359 


130±8.48 






HOY 


J053515 


26 


052109 


6 


83 


813 


-5 


352 


117±6.73 






HOY 


J053512 


71 


052139 


5 


83 


802 


-5 


360 


216±12.9 






HOY 


J053531 


30 


053308 


9 


83 


880 


-5 


552 


0.92±0.02 




11 


HOY 


J053521 


65 


053447 


1 


83 


840 


-5 


579 


0.38±0.01 




11 


HOY 


J053506 


50 


053335 





83 


776 


-5 


559 


5.10±0.17 




1/0 


HOY 


J053508 


52 


053558 


1 


83 


785 


-5 


599 


2.62±0.05 




1/0 


HOY 


J053504 


19 


053512 


1 


83 


768 


-5 


587 


1.37±0.11 




1/0 


HOY 


J053440 


92 


053144 


8 


83 


669 


-5 


529 


7.70±0.05 


13.3±0.27 


1/0 



a The two symbols ;:; and ... indicate that a source is extended or not detected, 
respectively. 

b Spectral energy distribution class following the scheme of I Greene et al. I il994l '). 
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Table 2. High-reliability flux time-series measured at 70 fim and 160 /mi between February and 

April 2011. 



Source ID 






Flux (Jy) 






A Flux a 




Feb. 26 


Mar. 6 


Mar. 14 


Mar. 21 


Mar. 31 


Apr. 10 


(%) 




PACS 70 


HOY J053529.44-045851.6 


6.36±0.07 


6.21±0.07 


6.20±0.05 


6.58±0.06 


6.35±0.05 


6.29±0.05 


6 


HOY J053530.33-045938.6 


0.96±0.01 


0.90±0.01 


0.85±0.01 


0.96±0.01 


0.98±0.01 


0.91±0.01 


15 


HOY J053529.50-045952.7 


12.8±0.10 


13.4±0.09 


12.2±0.09 


12.6±0.10 


13.2±0.09 


13.2±0.09 


10 


HOY J053515.98-050004.7 


7.34±0.07 


7.44±0.07 


7.20±0.06 


7.06±0.07 


7.19±0.08 


7.50±0.08 


5 


HOY J053518.12-050035.6 


27.3±0.17 


27.5±0.22 


27.5±0.19 


27.2±0.19 


28.2±0.21 


26.4±0.18 


7 


HOY J053522.27-050116.8 


31.0±0.19 


31.3±0.22 


32.3±0.21 


35.3±0.25 


37.0±0.28 


34.2±0.22 


19 


HOY J053523.26-050132.6 


51.8±0.32 


53.5±0.36 


49.6±0.32 


50.5±0.38 


55.1±0.38 


54.7±0.33 


12 


HOY J053519.61-050453.7 


3.58±0.04 


3.48±0.04 


3.49±0.04 


3.58±0.04 


3.77±0.06 


3.70±0.05 


8 


HOY J053525.89-050544.2 


11.2±0.09 


12.4±0.11 


12.5±0.10 


11.7±0.11 


12.8±0.10 


13.0±0.10 


17 


HOY J053525.51-050759.1 


10.5±0.07 


10.6±0.08 


10.5±0.07 


10.0±0.09 


10.5±0.07 


10.7±0.09 


7 


HOY J053524.23-050831.9 


7.61±0.08 


7.93±0.09 


7.25±0.10 


6.53±0.10 


7.04±0.09 


7.05±0.10 


21 


HOY J053524.71-051030.6 


93.6±0.90 


93.4±1.00 


95.5±1.13 


94.9±1.21 


97.7±1.15 


94.7±1.15 


7 


HOY J053527.56-050934.5 


460±6.08 


459±5.98 


457±5.99 


460±5.48 


431±6.51 


432±5.86 


8 


HOY J053523.29-051203.0 


50.7±0.44 


49.5±0.57 


46.2±0.39 


46.1±0.38 


50.7±0.48 


47.0±0.52 


11 


HOY J053520.15-051317.4 


45.2±0.34 


46.7±0.32 


46.2±0.30 


47.3±0.35 


48.1±0.30 


47.0±0.38 


7 


HOY J053519.47-051534.7 


27.0±0.57 


27.0±0.59 


26.4±0.55 


28.2±0.59 


27.1±0.56 


27.0±0.52 


7 


HOY J053515.69-052040.6 


191±6.43 


181±6.60 


190±6.27 


190±5.83 


204±6.45 


180±6.45 


17 


PACS 160 /im 


HOY J053529.50-045952.7 


20.9±0.68 


22.6±0.66 


21.7±0.63 


21.2±0.66 


21.4±0.69 


21.9±0.64 


7 


HOY J053518.12-050035.6 


37.3±1.54 


38.3±1.22 


39.2±1.45 


39.6±1.39 


39.4±1.60 


38.3±1.41 


11 


HOY J053522.27-050116.8 


67.4±1.34 


67.9±1.29 


71.5±1.47 


70.3±1.42 


72.7±1.63 


72.8±1.50 


12 


HHOY J053523.26-050132.6 


151±1.97 


156±2.07 


159±2.20 


160±2.25 


158±2.58 


166±2.13 


4 


HOY J053525.89-050544.2 


38.6±0.97 


39.5±0.64 


42.2±0.95 


42.4±1.01 


43.8±1.02 


42.9±1.0 


13 


HOY J053527.56-050934.5 


269±9.30 


279±9.35 


281±9.42 


282±9.28 


289±8.63 


274±8.69 


11 



a Peak-to-peak flux variations. Italic values indicate sources that pass our variability criterion (see Section [3] for 
details). 
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to be also variable in the mid-IR (left column of figures [2] and [3]) . However the shape of Herschel 
and Spitzer light curves are not directly comparable since they are not contemporaneous^] and they 
have significantly different observing cadences. 

The remaining 9 sources that have reliable light curves show flux variations in the 5 - 10% 
range. Although those sources exhibit variations below the 10% threshold, we cannot exclude 
that they might as well be intrinsically variable in the far-IR. In fact, based on the Herschel 
dataset currently available, we are only sensitive to variations with timescales in the range 10 - 
50 days, and we miss all variations on shorter and longer timescales. For instance, the light curve 
of HOY J053524. 71-051030. 6 does not satisfy our variability criterion according to its 70 /im flux 
only, however its Spitzer light curve shows clear variability in the mid-IR with a period of ~4 days 
(bottom right plot of figure [3]). If the mid- and far-IR variability had the same origin, and thus 
showed similar timescales at both wavelengths, then the sampling frequency of 10 days would not 
be appropriate to detect such rapid variations in Herschel bands. 



Although YSOs have been monitored on daily to yearly timescales from the optical t o the 
mid-IR and in all evolutionary stages (e.g. iHerbst et aDhoQcl : iMorales-Calderon et alJboilh . the 
variability parameter-space is still sparsely sampled in the far-IR. This is particularly true for high 
cadence monitoring and embedded protostars, as most far-IR m onitoring programs have targeted 



outbursting YSOs with fading timescales of several years (e.g. iKospal et al.1 120071 ). The present 



Herschel time-series observations thus fill in this parameter-space gap, and they show that the 



far-IR emission, which is a good tracer of the internal luminosity of protostars (|Dunham et al 



2008), can vary noticeably on timescales as short as a couple of weeks. This is orders of magnitude 
shorter than the dynamical timescales of the far-IR emitting material that orbits at hundreds to 
thousands of AUs from the central protostar. Such short timescales indicate that the mechanism(s) 
responsible for the observed variability occurs on smaller spatial scales, either at the surface or 
within the disk, or close to the central protostar (r<l AU) where timescales are consistent with 
weekly to monthly timescales. 



Non-steady accretion likely plays a role in the observed variability. IZhu et al .1 (|2009l . and refer- 
ences therein) argue that the typical accretion luminosity observed in protostars implies an accretion 
rate considerably lower than the predicted time-average infall rate in YSOs. This 'luminosity prob- 
lem' can be explained if infalling material accretes sporadically causing major accretion events that 
ar e sufficiently short- l ived s o that protostars are usually observed in quiescence. Numerical models 
of IVorobvov h Basul (|2010l ) manage to reproduce these accretion events for Class I and Class II 
YSOs when disk fragments fall onto the central protostar. Such energetic outbursts have been 
obser ved in FU-Orionis-type objects with typical timescales of 1 - 10 years (jHartmann Kenyon 
1990 ). Similar events might occur on shorter temporal and less energetic scales. 



3 There was no significant overlap in the Spitzer and Herschel visibility windows of Orion in 2011. 
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For instance, Flaherty et al.l (|201ll ) measured accretion rates from Pa/3 and Br7 observations, 
and they find that the accretion rate of the Class II object LRLL 31 varies by a factor of five, 
with the largest changes occurring on a weekly timescale. These episodic accretion events could be 
pictured as knotty filaments of gas funneling fro m the inner regio n of the circumstellar disk to the 
central star, likely through magnetic field lines (|Shu et al.lll994l ). and releasing packets of energy 
as the clumps hit the surface of the protostar. For embedded protostars, the energetic photons 
generated in the accretion shock are reprocessed by the dense and dusty surrounding envelop, and 
they eventually escape the protostar in the far-IR regime. The accretion luminosity variations that 
originate from episodic accretion events in the central region of a protostar propagate outwards 
and temporarily warm up the envelop, which in turn leads to detectable far-IR flux variations. 
The typical timescale for photons to reach the far-IR emitting material in the envelop (located at 
~ 10 2 — 10 3 AU from the central star) is roughly the light travel time, which is a few days, similar 
to the typical timescales observed in the light curves of Figure [2j 



The alternative scenario proposed by iFlaherty &; Muzerolld ()2010l ) might also cause a mea- 
surable far-IR variability in YSOs. It consists in a variable scale height of the disk inner edge 
that casts a shadow on the outer disk, thus cooling down the mid- and far-IR-emitting disk. This 
non-axisymetric disk model was used to reproduce the mid-IR variability observed in the YSO 
LRLL 31, and i t can potentially give constraints on the inner disk structure. In a subsequent work, 
Flaherty et al.l (|201ll ) have tested various models that may lead to variable scale heights of the 
inner edge of the disk, with origins ranging from variable accretion, perturbations by a companion, 
winds, and the influence of magnetic fields. Although LRLL 31 is more evolved than most objects 
observed in the present study, a similar mechanism could well explain the variations of the Class II 
source HOY J053530. 33-045938. 6 (bottom left plot in figure [2]). However, contemporaneous near- 
IR observations would be necessary to confirm the expected anti-correlated flux variations with 
respect to the far-IR light curves. 



The first 6 epochs of our monitoring program have demonstrated that protostellar emission can 
vary on relatively short timescales in the far-IR. Furthermore we find that the fraction of variable 
protostars in our sample is relatively high (>40%), though difficult to ascertain based on sparse 
70 fim light curves only. For a better sampling of the light curves, we have requested a higher 
observing cadence for the second Orion visibility window, and we have extended the initial span 
of the monitoring program with the last observations scheduled for autumn 2012, thus covering a 
period of nearly 2 years. These additional observations, with refined photometric measurements 
and a detailed SED modeling, should increase the number of variability detections and help us 
understand the nature of these objects and the mechanism responsible for the observed variability. 
A spectroscopic follow-up with Herschel/PACS is also scheduled for 2012 to monitor accretion ac- 
tivity in a couple of variable sources with atomic and mol ecular lines, particularly [01] , [CII], CO 



and H2O, already observed as part of the HOPS program (IMegeath et al.ll201ll : iPuravankara et al 

boul ). 
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This work is based on observations made with Herschel, a European Space Agency Cornerstone 
Mission with significant participation by NASA. Support for this work was provided by NASA 
through an award issued by JPL/Caltech. This research has made use of NASA's Astrophysics 
Data System. 

Facilities: Herschel Space Observatory (PACS). 
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Fig. 1. — Maps of the Orion A molecular ridge in the PACS 70 //in (top) and 160 p (bottom) bands 
displayed on a logarithmic scale. The first six visits to the targeted field were combined to create 
these artifact-free maps. Symbols on the top image indicate the location of detected protostars from 
table [TJ squares are Class II, circles Class I, and diamonds unclassified by the YSOVAR program. 
Thick symbols designate protostars with reliable light curves from table [2j Symbols on the bottom 
image show variable protostars following the definition given in section [31 
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Fig. 2. — Sample of reliable PACS light curves drawn from Table [2j The left column presents the 
light curves of variable protostars, and the right column contains those that show flux variations 
within the estimated photometric uncertainties (see section l2~3l for details). The set of graphs in the 
top row show the light curves at 70 /im and 160 /im when both are deemed reliable, while the other 
plots give 70 /im fluxes only. The horizontal dashed and dotted lines give the average flux of the 
sources and the ±5% variations around the mean, respectively, indicating our level of confidence 
for variability detections. The first epoch was obtained on Feb. 26, 2011, corresponding to the 
Herschel Operational Day 653, or MJD 55618. 



-14- 



HOY J053522. 27-0501 16.8 
IS0Y_J053522. 44-0501 14.3 



1.9 r 
1.8 h 

1.7 \ 
1.6 ' 

1.5 E* 4 

1.4 h 



10 20 

Time [MJD-551271 



HOY J053527. 56-050934. 5 
ISO Y_J05352 7.62-050932.9 



860 - 



840 - 
820 - 
800 * a 
780 - 



^ a^-^-.-a ^ 



760 - 



10 20 30 

Time [MJD-55127] 



HOY J053524.23-050831.9 
IS0Y-J053524. 30-050830. 8 



E 40- 



aA a^a a* 



20 

Time [MJD-55127] 



HOY J053515. 98-050004. 7 
ISO Y-J0535 1 6. 1 6-050002. 5 



r£ 115 - « a a 

E : « 



Time [MJD-55127] 



HOY J053530. 33-045938. 6 
IS0Y_J053530. 61 -045936.1 



HOY J053524. 71 -051030.6 
IS0Y_J053524. 73-051 030.1 




10 20 30 

Time [MJD-55127] 



Fig. 3. — Spitzer/IRAC light curves measured at 4.5 /xm on the 6 objects presented in figured! 
These observations were obtained as part of the YSOVAR program in 2009, they are not contem- 
poraneous with the Herschel data. The plots are arranged as in figure [2] with the YSOVAR name 
given in the plot title, and the horizontal dashed and dotted lines also represent the mean and ±5% 
variations, respectively. 



